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Highlights of earlier results by the author and others using the author’s Inside- 
Out. WAKE code on wake structures of LEO spacecraft are reviewed. For conducting 
bodies of radius large compared with the Debye length (large inverse Debye number), a 
high-Mach-number wake develops a negative potential well. Quasi neutrality is violated 
in the very near wake region, and the wake is relatively "empty" for a distance down- 
stream of about one-half of a "Mach number" of radii. There is also a suggestion of a 
core of high density along the axis. We report recent work on very large bodies in LEO. 

A comparison of rigorous numerical solutions with in-situ wake data from the AE-C 
satellite suggests that the so-called "neutral approximation" for ions (straight-line 
trajectories, independent of fields) may be a reasonable approximation except near the 
center of the near wake- This approximation is adopted here for very large bodies. 

In an earlier investigation of differential charging of small nonconducting 
bodies due to plasma flows, it was found that the scale of the voltage difference 
between the upstream and downstream surfaces ("front" and— I'wake" surfaces of a non- 
conducting body) due to a high-Mach-number plasma flow is governed by the ion drift 
energy. Hence kilovolt potential differences may occur in the solar wind, for example, 
between a Spacecraft and a piece of insulated material in its near wake. 

Recent work has concerned the "wake-noint" ootential of very large nonconducting 
bodies Such as the Shuttle Orbiter.. Using a cylindrical model for bodies of this size 
or larger in LEO (body radius up to 10^ Debye lengths), approximate solutions are pre- 
sented based on the neutral approximation (but with rigorous trajectory calculations 
for surface current balance). There is a negative potential well if the body is con- 
ducting, and no well if the body is nonconducting. In the latter case the wake sur- 
face itself becomes highly negative. The wake-point potential is governed by the ion 
drift energy. 


LARGE-BODY WAKE STRUCTURE: CONDUCTING BODIES 

Parker's wake-theory computer model for pillbox shapes (Inside-Out Method for 
warm ions - see refs. 1-3) was applied by the author and others in a number of wake 
calculations. High-voltage sheaths and wakes of large bodies require special numeri- 
cal techniques (see refs. 3 and 12 for generalization to 3-D geometries> CLEPH code). 


Wake of Moderately-Large Conducting Body in LEO 

First we present highlights of earlier results obtained (1976, see refs. 1-2) in 
a problem involving the wake of a large body in LEO, 100 Debye lengths in radius. The 
body is in the form of a disk oriented normal to the flow. For two cases (figs, la 
and lb) the parameter values are: 
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Despite possible inaccuracies, the consistency is such that physical conclusions 
may be drawn as follows.- In- this case the wake is seen to remain empty further down- 
stream than in the M=4 case. In- addition-* the suggestion is much- stronger that there 
is a central core of ambient density for both- ions and-.electrons along the axis. 
Moreover, the potential well is wider and longer than in the M=4 case* although the 
depth is about the same.- The normalized wake-surface fluxes are.? .-4 x 10 “30 (»'&") 
and. 4.2 x 10 ” 3 ° (final) for ions-, and 3-7 x 10“3 for electrons. The electron flux is 
slightly less than the M=4 value ,_and is again less than exp(-4). 

The conical region behind the disk where quasineutrality breaks down is now 
longer than in the M=4 case , extending to between zs4 and z=5 radii along., the axis . 

The region of wake disturbance is probably longer than 6 radii downstream, as in 
the M=4 case, but may not extend beyond about 2 radii in the transverse direction. 


Theory-Experiment Comparison for Ai'-C Satellite 

Next, we note that Parker's wake theory computer model has been applied by Samir 
and Fontheim (ref. 4) in a comparative study of ion and electron distributions in the 
wakes of ionospheric satellites. From a comparison between the theory and ion mea- 
surements on. the AE-C satellite, Samir and Fontheim show that theory and experiment 
agree fairly well in the "angle-of-attack" range between 90° and 135° • (The upstream 
and downstream directions are defined by 0° and 180°, respectively.) A significant 
finding is the fact that in that angular range even the "neutral approximation" for 
ions (straight-line trajectories, independent of electric fields) gives fair agree- 
ment with the measurements, (in the near-wake maximum rarefaction zone near 180 , 
both the neutral approximation and the self-consistent solution underestimate the 
measured ion densities - inferred from probe currents - by orders of magnitude i 
Electron data obtained by the Explorer 31 satellite also shows an underestimation 
near 180° by the Parker wake theory, although less pronounced.) 

The largest ratio of body-radius-to-Debye-length (that is, the inverse of the 
Debye number) treated by Samir and Fontheim (ref. 4) is Rj)=l62, in one of the AE-C 
cases. 

Figures 2a, b (from ref. 4) illustrate the geometry of the AE-C ion measurement, 
and the ion results for inverse Debye number l62. The locations of the ion current 
observation points, and of the numerical grid points at which densities were calcu- 
lated, are shown in figure 2a. The geometry of the theoretical model is that of a 
pillbox cylinder with its axis parallel to the flow, while the true geometry is that 
of a pillbox cylinder in a "cross-flow," that is, with its axis perpendicular to the 
flow. In spite of this, the theory-experiment comparison is deemed by Samir and 
Fontheim to be meaningful, in view of uncertainties in the calculations and estimated 
measurement errors. (The depth in the direction of the flow is the same for both the 
satellite and the model, and the cross sections presented to the flow are nearly the 
same.) The current probe moves on a circular arc at a radial distance of about 1.5 
satellite radii. 

In figure 2b, the measured angular profile is shown together with the neutral 
approximation (zero-th iteration) and the self-consistent solution (15-th iteration). 
The self-consistent solution is closer to the experimental profile, in the angular 
range 90 ° - l4?°, than the neutral approximation. Near 180°, the self-consistent 
solution is 2 to 3 orders of magnitude below the measured data, while the neutral 
approximation is about 10 orders of magnitude lower. 


the F r heim •*«*• «-» 

the Inverse Debye number (ratio of body radius to Debve len^h^H m ° re f &CUrately as 

ma 

the quasineutrality regime, at sufficiency larg e Sref ^^Sre" 1 ” 8 ' 1 ^ 

Wake of Very Large Conducting Body in LEO: Re^nt Results 

previously^ ^ ^ any treated 

relatively expensive. However a r,»»«nn!hi * calcu . Latlon becomes computationally 

of the "neutral approximation" ’for ions. Tha^is^the^n 0 + S a f f ° rd ? d through the use 
space charge densitv arp «« ^ iu- * ^ lon tra J^ctories governing ion 

fieldT Thf~ electron spLe charl de^i^ / Wera and unaffected by the 

factor", that is, the exponential tTtht vll, as ? ujned . to b ? S iven b 7 the "Boltzmann 
extent this appr^ximat^X ^ potential. To some 

discussed above. In anv case it rtl , a1 f + +• , lr and Fontheim m-situ comparison 

sights with a minimum of computational expense ^ifa^rox^ t0 physical in ‘ 

et al (ref. u,. (We compute current tal^ate^' ^“t^eSi* f* 

a ion?XS“^ns rl ^^ th tVS not s f ellIte - ln or 

potential of 3 kT/e, is shown in fieur-M - 3 n v . assumed to have a dimensionless 

numbers ranging from 10 to 10*5, and flow Mach nSers 2 % ^8 Wi ** inve f e » eby e 
the wake potential profile varies with i-mmro nambers 2 ’ 5 and 8. Eigure 3a shows how 
8. The profiles for inveSSbye Sere W '"•Rf” 4 h n ™ ter * 

obtained earlier for a sphere hv Kji^i ’ „. and are Slmila r to results 

(ref. 11) results are for inverse Debye n^berfup'to 1<? 'J''. 11 ’’ ?". Klel St al 
tions to lO^. The wake DotervHAi u ^ have extended the solu?^- 

numbers greater than about 10.. The magnitude of’thi^'” 1 ” 1 ^™ f° r lnverse Deb J' e 
19, respectively, for inverse Debye Sere 10* m3 10 ' U and 

how the wake potential profile varies with Marh’ and 10 . Figure 3b shows 

ber = io5>. The depth of the SJIS? V * • number, for fixed inverse Debye num- 

Mach number and inverse Debye numbed Pi'gmTso^hovs «Tf b ° th i ’ lcrt ‘ a1 "* 

number = 8 and. inverse Debye number = iof. 3 h e( I ul Potential contours for Mach 

about^P! rrnS7coS?“rbo^° £i!5 t, s^s t s. (ta T" r ye Mb t r 

covered with nonconducting tiles Hence it milt hJ + f , surface (about 97%) is 


WAKE STRUCTURES AMD DIFFERENTIAL CHARGING OF SMALL AND LARGE 
NONCONDUCTING BODIES DUE TO PLASMA flows 


Differential Charging 

partly or entirely insSIting C and S the charee^^^^? 11 thS spacecraft surface is 
Point over the surface. In ^aSlit 


238 


insulated area, due to electrons escaping from it the sunlit area tends to become 
positively charged relative to the surrouhding dark areas (refs. 5-7). Another mech- 
anism of differential charging, which is less familiar and appears to have been 
treated only very, recently (ref. 8), is that due to the relative motion between a 
nonconducting spacecraft and the external plasma (e-.g., a spacecraft in- the iono- 
sphere or in the solar wind). The fluxes of ambient ions and electrons on- the wake 
surface are not the same as on the front surface. For high velocities of relative 
motion compared with the mean ion thermal velocity, whether this occurs in the iono- 
sphere (due principally to spacecraft motion) or in the solar wind (due principally 
to plasma motion), there is a significant differential in the ion fluxes, but a neg- 
ligible differential for the electrons.. Since the net current density must vanish 
locally at each surface point in the steady state, this plasma- flow effect leads to a 
larger negative equilibrium potential on the wake Surface than on the front surface.. 
If there is photoemission as well on the front surface (as in the solar wind), this 
differential charging is enhanced . As shown below, this plasma- flow effect can gener- 
ate differences between, the front and wake surface potentials amounting to many kT/e 
(where T is the temperature, k is Boltzmann’s constant, and e is the electron charge), 
together with a potential barrier for electrons. The potential difference can be 
expected to be of the order of volts in the ionosphere, and one kilovolt in the solar 
wind, that is, of the order of the ion drift energy (ref. 8). 

Even weak, differential charging can interfere with measurements of, say, weak 
ambient electric fields or low-energy particle spectra, and it can create electron 
potential barriers which can return emitted photoelectrons or secondary electrons to 
the surface and lead to erroneous interpretations of the data (ref. 9). This type of 
electron potential barrier is distinct from, and- should not be confused with, the 
more familiar space-charge potential minimum which can be produced by emitted- 
electron space charge (ref.. 10) and is not due to differential charging. The barrier 
produced- by differential- charging effects may be mo re import an t than the potential 
minimum caused by space charge. 

The next section results show what may be expected: (a) in the ionosphere for 

small insulated objects, small meteroids, or small parts of a spacecraft (e.g., a 
painted antenna) located within the wake region of a moving spacecraft, and (b) in 
the solar wind for an entire spacecraft, or small natural bodies in the solar system. 
Following the next section, the wake structure and differential charging of very 
large nonconducting bodies in Low Earth Orbit will be treated. 

Differential Charging of Small Nonconducting Body 

In the problem treated next (see fig. 4), we assume the nonconducting spacecraft 
to have a "pillbox" shape, and to be in a flowing plasma, with the plasma flow along 
the axis, from the "front" region toward the "wake" region. The plasma is taken to 
be ionized hydrogen and is assumed to have a velocity of flow 4 times larger than the 
most probable ion thermal velocity (ion "Mach number" =4). (in the solar wind, this 
Mach number would be approximately 10.) Since the unperturbed ion flux to the wake 
surface is about 9 orders of magnitude smaller than the corresponding ion flux to the 
front surface, and since the electron fluxes are about the same to the front and wake 
surfaces, there will be a significant differential between the equilibrium potentials 
at the front and wake surfaces (see below). 

Using the Inside-Out Method, current densities of ions and electrons are evalu- 
ated at many points on the spacecraft surface (refs. 7-8). The local surface poten- 
tials were varied until current balance was achieved at each point. 
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Wake Structure and. Differential Charging of Very Large 
Nonconducting Bodies- in- LEO Plasma Plows : Recent-Results 

There is considerable interest in the charging and electric fields of the Shuttle 
Orbiter. This is an important example or a very large spacecraft in Low Barth Orbit 
(Inverse Debye- number about loM with most of its Surface (about 97%) nonconducting . 
Only the small area in the vicinity of. the engines Is conducting and electrically 
grounded to the main frame. Figures 5a and 5b indicate how the Orbiter may be sub- 
jected to different types of differential charging depending on. its orientation with 
respect to the plasma-flow direction. In figure 5a, the Orbiter- is moving "nose- 
forward," i.e.. heading into the flow. The wake-point potential (location indicated 
by a cross) occurs essentially in the engine area, and thus defines also the Orbiter's 
ground potential. The rest of the spacecraft surface is electrically isolated and 
has in general a different potential distribution. The cargo bay area is a "side" 
region according to the terminology of the previous section. In figure 5b the 
Orbiter is moving "belly- forward." With this orientation the wake-point potential 
occurs in the cargo bay area ^ which is electrically isolated from the Orbiter ground. 
The ground is defined by a different potential attained by the engine area. In the 
shown- orientation, the engine area is a "side" region. 

Hence, the maximum negative ground potential of the Orbiter would occur when the 
Orbiter is in the nose-forward orientation, while the cargo-bay potential would be 
intermediate between this and the plasma potential. With the belly-forward orienta- 
tion, the roles of ground potential and cargo-bay potential would be reversed, with 
the cargo bay at maximum negative potential, and Orbiter ground at intermediate 
potential . 

Iil the present paper the wake structure and the wake-point potential of a very 
large nonconducting body in LEO such as the Orbiter are calculated- using certain 
approximations- The geometry is modeled by a circular cylinder as illustrated in 
figure 6, The wake point is the isolated area indicated by a cross in the figure. 
Again, because of computational expense, we use the neutral approximation, hut only 
for ion space charge. However, the differential charging, e.g. the wake-point poten- 
tial, is calculated rigorously by current balance using Inside-Out-Method trajectories 
(refs. 7-8), for both ions and electrons, in the resulting electric field distribution. 

For a nonconducting body of any size, current balance at the wake point results 
in significant negative wake-surface potentials. (Nonconducting bodies were not 
treated by Kiel et al.) Figures 7a and 7b show results for inverse Debye number 10^, 
and Mach numbers 2, 5 and 8. There is no potential minimum. Instead the wake point 
attains the highest negative potential, resulting in a monotonic rather than non- 
monotonic potential profile in the wake. Figure 7a shows how the wake-point poten- 
tial increases with increasing Mach number, for a fixed inverse Debye number = 10^. 

The wake-point potential magnitude is about 8, 20 and 36 kT/e, respectively, for Mach 
numbers 2, 5 and 8. Figure 7b shows equipotential contours for Mach number = 8 and 
inverse Debye number = 105. These contours (nonconducting body) may be compared with 
those of a large conducting body with the same parameters (fig. 3c). 

fable 1 shows how the wake surface potential of a nonconducting large body 
varies with Mach number and inverse Debye number. Evidently, the wake, surface poten- 
tial is insensitive to inverse Debye number. The table also gives the values of the 
dimensionless current density (equal of course for ions and electrons) at the wake 
surface. For comparison, also shown are the ion currents that would result from 
using the neutral approximation to calculate currents (see ref. 7). These are seen 


to be many orders of magnitude smaller than the more realistic currents calculated 
using ion trajectories affected by the field. 

For large nonconducting bodies in hi gh-Mach-number flows, the wake-to-front 
potential difference generated by the flow is less than but of the order of the 
potential-equivalent of the ion drift energy.- This result is similar to that 
obtained above for the case of a small nonconducting body. 

Finally, we illustrate in figures 8a, 8b and 8c examples of intricate 3-D large- 
body geometries of aerospace interest (including the Orbiter) for which a wake- 
modeling capability will be achieved using techniques presently under development at 
Lee W. Parker, Inc-. 
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Figure 1. - Concluded. 












SATELLITE 


(a) Theoretical model (solid) versus real satellite geometry (dotted). (Dots 
denote ion current observation points; X's denote numerical arid points 
at which densities are calculated.) 
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(b) Measured angular profile on AE-C satellite (large body; 162 Debye lengths) 
conroered with neutral-approximation theory (iteration zero) and self- 
consistent theory 0 eration 15). 

Figure 2. - Geometry of AE-C ion measurements. 
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(a) Variation with inverse Debye number (b) Variation with Mach number, at fixed 
at fixed Mach number *8. inverse- Debye number ■ 10 . 



(c) Equipotential contours. Mach number = 8. Inverse Debye 
number ■ 10^. (The point marked "x" is the position of 
peak potential « 19; dimensions in units of spacecraft radius.) 

'igure 3. - Wake potential profiles (dimensionless potential) ar.d equi- 
potential contours in wake of conducting cylinder with 3 kT/e surface 
potential. $ = potential in units of kT/e; r = downstream distance in 
units of spacecraft radius. 
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WAKE BOUNDARY 




POTENTIAL 


(a) Nose-forward orientation. 


(b) Belly-forward orientation. 


F Arbiter aroInd^M^j^'' *? EE0 p ' asnB Indicating wake points and 
g ound potential points. Very large nonconducting body. 
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(a) Variation with Mach-mwnber at fixed inverse Debye number of 10^. 



(b) E 5“jP°tentia1 contours (dimensionless potential) in wake of nonconducting 
cylinder. Surface potential distribution from 3 to 36.2, in units of 

nlh»i det ® rmined Pointwise current balance. Mach number = 8. Inverse 
Debye number > 10 . (Dimensions in units of spacecraft radius). 

Fl ?infn.^c” Wake b POt f ntial P rofiles (dimensionless potential and equipotential 
contours in wake of nonconducting cylinder. * . potential in units of 
Ki/e; r * downstream distance in units of spacecraft radius. 
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(b) C-130 Hercules aircraft model. 

Fl patches.~ Three " dlmensional computer models constructed of quadrilateral 
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(c) Starship Enterprise model 
Figure 8. - Concluded. 


252 




Art-' 


,v ~y 


• &». 



